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Abstract 

An extensive study of GeV light-ion-induced multifragmentation and its possible 
interpretation in terms of a nuclear liquid-gas phase transition has been performed 
with the Indiana Silicon Sphere (ISiS) 47r detector array. Measurements were per- 
formed with 5-15 GeV/c p, p, and vr" beams incident on ^^^Au and 2-5 GeV ^He 
incident on ""*Ag and ^^^Au targets. Both the reaction dynamics and the subse- 
quent decay of the heavy residues have been explored. The data provide evidence 
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for a dramatic change in the reaction observables near an excitation energy of E*/A 
= 4-5 MeV per residue nucleon. In this region, fragment multiplicities and energy 
spectra indicate emission from an expanded/dilute source on a very short time scale 
(20-50 fm/c). These properties, along with caloric curve and scaling-law behavior, 
yield a pattern that is consistent with a nuclear liquid-gas phase transition. 

Key words: Att detector array,multifragmentation, reaction dynamics, nuclear 
density, reaction time scale, caloric curve, scaling laws, liquid-gas phase transition 
PACS: 25.70.Pq, 25.55.-e 
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1 Introduction 



One of the most important signals of the formation of hot nuclear matter is 
the emission of nuclear clusters, or intermediate-mass fragments (IMF: 3 < Z 
< 20). From studies of the IMF yields in energetic nuclear reactions one hopes 
to gain greater insight into the thermodynamics of highly-excited nuclei and 
the nuclear equation-of-state at low densities. An important aspect of such 
studies is the identification of a possible nuclear liquid-gas phase transition 
[1,2,3]. 

IMF emission was first observed in the 1950 's, when beams of protons and 
alpha particles with energies in the GeV range became available. Emulsion 
and radiochemical measurements of the reaction products showed that the 
probability for cluster emission increased strongly with beam energy, sug- 
gesting their association with the decay of highly excited nuclei. The emul- 
sion measurements also provided evidence for multiple fragment emission, or 
multifragmentation. These data led to the concept that the reaction mecha- 
nism could be schematically viewed in terms of a two-step intranuclear cas- 
cade/statistical emission model, with IMFs emitted primarily in the evapora- 
tion stage. One complication with this model was the observation that lighter 
clusters exhibited forward-peaked angular distributions, suggesting that a 
prompt non-statistical mechanism must also be present. Another perceptive 
proposal was that pion production and reabsorbtion could be a major mecha- 
nism for excitation-energy deposition in the heavy residual nucleus[4,5]. These 
studies, carried out at Lawrence Berkeley Laboratory, Brookhaven National 
Laboratory and in Russia, laid the groundwork for future studies and are 
reviewed in [6,7,8,9]. 

The development of silicon semiconductor detectors in the 1960's made it 
possible to perform measurements of inclusive IMF yields and kinetic energy 
spectra. Bombardments of heavy targets with 5 GeV protons provided a more 
systematic understanding of the earlier work, and showed clearly the exis- 



tence of two mechanisms for IMF production - one equilibrium-like and the 
other a fast nonequilibrium process [10,11]. An important aspect of the spec- 
tra was a downward shift in the apparent Coulomb barrier relative to lower 
energy reactions. It was suggested that this shift was due to a modified density 
distribution, perhaps due to an expanded source. Later, similar inclusive mea- 
surements with protons up to 350 GeV [12] confirmed radiochemical studies 
that indicated the IMF emission probability reaches a maximum near 10 GeV 
and remains constant thereafter [13]. These later measurements also tracked 
the evolution of the Coulomb peak displacement with beam energy, and found 
its onset to be near 4 GeV. This result was interpreted as a possible sign of 
critical behavior [1,14], and stimulated widespread interest in this possibility. 

To search for more conclusive evidence of a phase transition and critical be- 
havior, it was necessary to develop large-solid angle detector arrays to provide 
fragment kinetic-energy spectra, multiplicity information, event topology and 
calorimetry. Among the most important features of such a device are: (1) nu- 
clide (Z and A) identification of all products, including neutrons, (2) spatial 
characterization with good granularity; (3) low detection thresholds; (4) good 
energy resolution over a large dynamic range and (5) efficient, reliable de- 
tector calibration techniques. Another important factor in detector design is 
the geometry imposed by reaction kinematics; e.g. IMF fragment distributions 
are nearly isotropic in light-ion-induced reactions, with a small component in 
the beam direction, whereas they are strongly forward-focused in heavy-ion 
reactions. 

With improvements in detector and data-acquisition techniques, it became 
possible to construct complex detector arrays for performing exclusive mea- 
surements that met most, but not all, of the above conditions for complex 
fragments. The first such array was the LBL Plastic Ball/Wall [15], which 
demonstrated the existence of events with high multiplicities of IMFs and 
light-charged particles (LCP: Z = 1,2) [16]. Later measurements at LNS Saclay 
confirmed the high IMF multiplicities and yield dependence on bombarding 
energy [17]. These measurements also showed that the IMF spectral peak en- 
ergies decreased as a function of increasing multiplicity, consistent with the 
expansion scenario. Subsequently, several in detector arrays were constructed 
for the study of light-ion induced multifragmentation: the EOS TPC [18], the 
Berlin Silicon/Neutron Ball [19], the Dubna FASA array [20] and the Indiana 
Silicon Sphere (ISiS) [21]. At the same time heavy-ion accelerator technology 
had advanced to the stage where it became possible to study multifragmenta- 
tion in A + A reactions, which lead to the construction of several additional 
detectors [22,23,24,25,26,27]. At the outset it should be stressed that there is 
general concordance among the results of all of these projects with those from 
ISiS, as can be found in [27]. In the following sections, we survey the salient 
results of the ISiS program with GeV light-ion projectiles, and address the 
question: is there evidence for a nuclear liquid-gas phase transition? 



2 The ISiS Experimental Program 



The Indiana Silicon Sphere project was initiated in the late 1980's with a spe- 
cific focus on GeV light-ion-induced reactions on heavy nuclei. Light ions bring 
a unique perspective to multifragmentation studies in that they emphasize the 
thermal properties of the disintegrating residue, with minimal distortions of 
the spectra due to rotational and compressional effects. Experimentally, GeV 
light ion beams form only one emitting source, which decays in a reference 
frame that is very close to that of the source frame (source velocities ~ 0.01c). 
In addition, the energy-deposition mechanism produces a continuous distribu- 
tion of excitation energies, permitting broad coverage of the excitation energy 
distribution in a single reaction. 

Based upon previous inclusive and semi-exclusive studies [17], the ISiS Air 
array was designed in a spherical geometry, with very low detector thresholds 
and a large dynamic range for LCPs [LCP = H and He] and IMFs [21]. Particle 
identification was based on triple telescopes composed of an ion-implanted 
silicon detector, preceded by a gas-ion chamber for low energy fragments and 
followed by a Csl (T£) crystal for the most energetic particles, primarily LCPs. 
An overall schematic diagram of ISiS is shown in Fig. 1, accompanied by a 
photo in Fig. 2. The array consists of 162 particle-identification telescopes, 90 
covering the angles 14°-86.4° and 72 spanning 93.6°-166°. The telescopes are 
arranged in eight rings, each of which is composed of 18 truncated ion-chamber 
housings. In the forward-most ring, the Si and Csl detectors were divided into 
two segments to increase granularity. Detector acceptance was 74% of An. 

A schematic drawing of the detector telescopes is shown in Fig. 3. The de- 
tectors are mounted in gold-plated, high- conductivity copper cans. The first 
element in each telescope (AE) is an axial- field gas-ionization chamber (GIC), 
operated at ~ 200 V and a pressure of 15-20 Torr of CsFg gas. These con- 
ditions permitted identification of fragments with energies as low as ~ 0.8 
MeV/nucleon. All detectors operate in a common gas volume in each hemi- 
sphere, with vacuum isolation provided by a ~ 120-150 /ig/cm^ stretched 
polypropylene window covered with a thin graphite coating. The ion cham- 
bers are followed by an ion-implanted silicon detector of thickness 500 /im, 
which is sufficient to stop E/A ~ 8 MeV LCPs and IMFs. The GlC-Si tele- 
scopes provided Z resolution for Z = 1 - 16 fragments in the energy interval 
from E/A ~ 0.8-8.0 MeV. Due to the GIC energy resolution, mass identifica- 
tion was not possible in the GlC/Si pair. The final element of each telescope 
was a 28 mm-thick CsI(T£) crystal, read out by a photodiode. These crystals 
provided an energy acceptance of 1 MeV < E/A < 92 MeV. The Si-CsI(Tl) 
pair provided Z and A identification for 8 MeV < E/A < 92 MeV particles 
and "grey particle" detection for fast protons and pions up to ~ 350 MeV. 
The ISiS telescopes did not detect neutrons or Z-identify heavy residues and 




Fig. 1. Assembly drawing of the ISiS detector array. Components are as follows: (1) 
center support ring; (2) gas-vacuum separation window; (3) arc support bars; (4) 
partition disks; (5) beamline support cones; (6) target ladder assembly; (7) steel 
rails for opening housing cans and (8) vacuum chamber. 



fission fragments. 

Detector signals were sliaped and amplified by means of charge-sensing pream- 
plifier/linear shaper NIM units, with gains custom-designed for each detector 
type. Analog signals were digitized by 12-bit 16-channel peak-sensing ADCs 
and fast signal discrimination and multiplicity sensing was accomplished with 
16-channel discriminators and time-to-digital converters. Voltages for each de- 
tector type were supplied by in-house-designed, computer-controlled bias sup- 
ply units. 



The event hardware trigger was generated from the fast outputs of the dis- 
criminators and required a minimum of three detected particles. This criterion 
was imposed by the high backgrounds associated with synchrotrom accelera- 
tors and biased the data against low-excitation energy events (E*A < 1 MeV). 
In software, only events with three thermal particles, one of which with Z > 2, 
were accepted. More complete details for the ISiS array can be found in [21]. 




Fig. 2. Photograph of the ISiS detector loaded with detector telescopes. 

Three inultifragmentation campaigns were carried out with ISiS: (1) E228 with 
1.8 -4.8 GeV ^He ions at the Laboratoire National Saturne in Saclay, France; 
(2) E900 at the Brookhaven AGS accelerator with 5.0-14.6 GeV/c proton and 
TT^ beams; and (3) E900a at AGS with tagged 8.0 GeV/c antiproton and vr^ 
beams. The number of events with each beam is summarized in Table 1. 

In the following sections we present the results obtained with ISiS. We first 
examine the collision dynamics, then the thermal observables, and finally the 
thermodynamics and scaling-law behavior of the data. 



3 Reaction Dynamics 



3.1 Excitation Energy Deposition 



In GeV light-ion-induced reactions the dissipation of radial beam energy into 
internal excitation of the target-like residue proceeds through a complex fast 
cascade of nucleon-nucleon collisions. This mechanism is abetted by the ex- 
citation of A and higher resonances, followed by the reabsorbtion of some of 
the decay pions [28,29,30,31,32]. For antiproton beams, excitation-energy de- 
position is further enhanced by the reabsorbtion of some fraction of the 4-5 




Fig. 3. Drawing of an ISiS arc bar segment for the forward hemisphere with the 
angular coverage of each telescope labeled. Each segment is part of an 18-member 
ring; the forward- most element is divided into two segments. Forward-angle rings 
are identified as follows: 14-22°; 22-33°; 33-52°; 52-69°; 69-86.4°. Backward-angle 
rings are: 93.6-111°; 111-128°; 128-147°; and 147-166°. 

annihilation pions [33]. Although the energy dissipation process is relatively 
inefficient, the cascade step is capable of imparting up to ~ 2 GeV of excitation 
energy in heavy nuclei. Equally important, randomized/equilibrated residues 
are produced with a continuous distribution of excitation energies, essentially 
providing nearly a complete excitation function in a single reaction. 

Fig. 4 depicts the predictions of a BUU (Boltzmann-Uehling-Uhlenbeck) cal- 
culation [34,35] that traces the time and density evolution of a central (b = 
2.0 fm) collision between a 14.6 GeV/c proton and a ^*^^Au nucleus. Initially, 
a local density depletion develops along the projectile trajectory as forward- 
focused nucleons are ejected on a fast time scale. After about 30 fm/c the 
nuclear matter density becomes more uniform, indicating a random distribu- 
tion of nucleons. However, the average density is observed to be significantly 
lower than that of the original target, creating a hot, dilute nucleus. As the 
reaction time increases, the angular distribution of the emitted particles be- 
comes more isotropic. At this point the distinction between low-energy cascade 
nucleons (nonequilibrium) and evaporative nucleons becomes blurred. Also, it 
is significant to notice that the heavy residue trajectory has a component 
transverse to the beam direction. 



Table 1 

Number of events analyzed for each system (in parentheses). Event acceptance re- 
quires at least three thermal charged particles in sihcon detectors, one of which 
must have Z > 2. 



Beam 


Target 
197Au 


5.0 GeV/c 


Energy/Momentum 


9.2 GeV/c 


8.0 GeV/c* 


8.2 GeV/c 






(1.0 X 10^) 


(2.5 X 10^) 


(2.4 X 10^) 


(1.4 X 10^) 


P 


197Au 


6.2 GeV/c 


9.2 GeV/c 


12.8 GeV/c 


14.6 GeV/c 






(2.4 X 10^) 


(1.7 X 10^ ) 


(1.4 X 10^) 


(1.1 X 10^) 


P 


197Au 




8.0 GeV/c* 
(5.5 X 10^) 






^He 


natp^ 


1.8 GeV 


3.6 GeV 


4.8 GeV 








(4.9 X 10^) 


(3.0 X 10^) 


(1.9 X 10^) 






197Au 


1.8 GeV 
(4.0 X 10^) 




4.8 GeV 
(2.9 X 10^) 




*tagged beam 













In Fig. 5 the effect of entrance-channel beam momentum on excitation energy 
per nucleon (E*/A), average density (< p>/po), entropy per nucleon (S/A) 
and residue mass (A) is explored as a function of time. Calculations are for 
an impact parameter b = 2.0 fm for the p + ^^^Au reaction at momenta 
of 6.2, 10.2, 12.8 and 14.6 GeV/c. It is observed that the excitation energy 
and entropy per nucleon increase with beam momenta, while the source mass 
and average density decrease. At long reaction times there is little difference 
in E*/A and < p>/po indicating a saturation in these variables. Also, in 
all cases the entropy per nucleon remains nearly constant beyond 30-40 fm/c, 
consistent with the existence of a randomized system. The most significant de- 
pendence on increasing beam momentum is the systematic decrease in residue 
mass, or in terms of the emitted particles, a greater contribution to the cas- 
cade/nonequilibrium yield. 

All p + A reaction dynamics calculations predict that the excitation energy 
probability decreases exponentially with increasing excitation energy. For this 
reason, it is of interest to investigate the relative efficiency of different hadron 
probes in depositing excitation energy. Fig. 6 shows results of an intranuclear 
cascade calculation [31,36,37] for the average excitation energy (E* > 50 MeV) 
as a function of beam momentum for proton, negative pion and antiproton 
beams. For the p and vr" cases there is little difference, since the cascades 
follow similar paths. On the other hand, the annihilation pions from the p 
interaction enhance <E*> significantly. The inset in Fig. 6 demonstrates that 
the excitation-energy distribution for antiprotons is expected to extend to 
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Fig. 4. Nuclear density within the reaction plane XZ around the position coordinate 
Y=0, where Z is along the beam axis at different times in the p + ^^^Au reaction 
at 14.6 GeV/c. Calculation is for b=2 fm, with the {ain-,U'^) transport simulations 
discussed in Sec. 3.2. [35] 

higher values than for protons and pions. 

In comparing the predictions of the dynamics codes with experimental data, 
the primary quantities of interest are the energy deposited in the statistical 
residue E* and its mass A. To do so, however, requires separation of cas- 
cade/preequilibrium emissions from those associated with the equilibrium-like 
heavy residue prior to event reconstruction. This calorimetry procedure is de- 
scribed in Sec. 4.1. 

In Fig. 7 we show the reconstructed probability distributions for excitation 
energy and residue mass for several systems studied in this work. The recon- 
structed E* distributions shown in the left panel of Fig. 7 demonstrate that 
the largest population of high excitation-energy events is achieved with the 8.0 
GeV/c p beam and the lowest with the 5.0 GeV vr" beam. Thus, the data are 
qualitatively in agreement with the INC calculations (intranuclear cascade) 
shown in Fig. 6, although the calculations extend to somewhat higher ener- 
gies than the data. The residue mass distributions in the right panel of Fig. 
7 show a different pattern. In this case the 14.6 GeV/c proton beam produces 
the lightest residues and the 5.0 GeV/c vr— beam the heaviest, a result rel- 
atively well reproduced by the calculations. This mass dependence on beam 
momentum can be understood as a consequence of the fast cascade, which 
produces an increasing number of fast knock-out particles as the beam mo- 
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Fig. 5. Excitation energy per nucleon of the source, mean density, entropy per 
nucleon and residue mass number, as a function of time in the p + Au reaction at 
b=2 fm and different incident momenta, from the {am, U^) transport simulations 
[35], discussed in Sec. 3.2. 

mentum increases [28,32]. This process produces the saturation in excitation 
energy observed for hadrons with momenta greater than ~ 8 GeV/c. That 
is, the increase in total beam energy available for E* deposition is counter- 
balanced by the the loss of energy due to mass loss AA during the fast cascade. 

The relative effectiveness of various beams in depositing high excitation ener- 
gies (Fig. 7) is emphasized in the bottom panel of Fig. 8. Included here are 
comparable data from the 4.8 GeV ^He + ^^^Au reaction [38] and from the 
1.2 GeV p + ^^'''Au reaction [39]. In order to emphazise the probability for 
forming highly excited systems, all data are normalized to probability P(E*) 
= 1 at E* = 400 MeV. 

Fig. 8 confirms that the 8.0 GeV/c p beam produces a significant enhancement 
of high excitation energy events, particularly in the multifragmentation region 
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Fig. 6. Intranuclear cascade predictions [31] of the average excitation energy for 
events with E* > 50 MeV are shown as a function of beam momentum for p, tt^ 
and p beams incident on ^^^Au. Inset compares the excitation energy probabihty 
distributions for 8 GeV/c vr^ and p beams 

above E* > 800 MeV. This figure further supports the predictions of the 
cascade code; i.e., enhanced E* with p beams, saturation above ~ 8 GeV/c 
beam momentum, and exponentially-decreasing probabilities for high E*/A 
values. The behavior of the ^He beam can be understood as due to its average 
beam momentum of 2.6 GeV/c per nucleon. 

Another perspective on the influence of reaction dynamics is provided by stud- 
ies with ^He beams [40,41]. Fig. 9 shows the probability distribution of the 
total observed thermalized energy and the correlation with the transverse frag- 
ment kinetic energies for the ^He -|- ""*Ag and ^^^Au reactions. The thermal 
energy, which is strongly correlated with excitation energy, is significantly 
lower for the lighter ""*Ag target. However, when source mass corrections are 
made, the energy per nucleon distributions are very similar for both targets 
[40]. The projectile-target effect on the saturation of excitation energy with 
beam energy is in good agreement with INC calculations [30,41] for this lower 
momentum, complex projectile. The lower right-hand frame of Fig. 9 shows 
that the probability distributions are identical for the ^He -|- ""*Ag at 3.6 and 



12 



10 



LjJ 

CL 



10 



10 




---- 8.0 GeV/c 71" 

8.0 GeV/c p 

• 5.0 GeV/c 7T" 
D 9.2 GeV/c TT" 
A 10.2GeV/cp 
T 14.6GeV/cp 




1000 

E* (MeV) 



100 



200 



A. 



Fig. 7. Measured excitation energy (left frame) and residue mass (right frame) 
probabilities for several of the systems studied in this work, as indicated on the 
figure. '^P(E*) = 1. Data for 6.2 and 12.8 GeV are not shown but are consistent. 
Values < 250 MeV are uncertain because of missing neutrons in the calorimetry 
procedure. 

4.8 GeV, indicating the onset of deposition energy saturation near 3.6 GeV. 
The slight difference at 3.6 and 4.8 GeV for the transverse energy correla- 
tions is due to the difficulty in removing nonequilibrium events from this sum, 
which illustrates the uncertainties involved with the use of transverse energy 
as a gauge of excitation energy depostion. 



3.2 BUU Simulations 



While the equilibrium-like events are of primary interest for multifragmen- 
tation studies, the nonequilibrium component of the spectrum is important 
for understanding the reaction dynamics. In order to investigate this facet of 
the data, predictions of a BUU code that includes d, t and ^He cluster for- 
mation have been compared with cascade/preequihbrium d/p, t/p and ^He/p 
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Fig. 8. Bottom: the probability for observing events with excitation energy greater 
than E* > 400 MeV relative to the probability for events with E* = 400 MeV. 
Systems are indicated on the figure. Top: average mass loss AA in the fast cascade 
as a function of excitation energy. Systems are defined in bottom frame. 

ratios for the p + ^^'^Au reaction between 6.2 and 14.6 GeV/c [34,35]. BUU 
calculations were performed for various reaction times with and without a 
momentum-dependent potential and with both free and in-medium cross- 
section options. Results for the 14.6 GeV/c p -|- ^^^Au reaction are shown 
in Fig. 10. From examination of the LCP ratios for all four bombarding ener- 
gies, the best agreement with all the data is found when both a momentum- 
dependent potential and in-medium cross-sections are employed in BUU code, 
with a most probable reaction time of t ~ 65 fm/c. From comparison with Fig. 
5, entropy considerations suggest that randomization of the nucleon momenta 
in the heavy residue occurs after about t ~ 30 fm/c, indicating an additional 
30-40 fm/c is required to reach a state of quasiequilibrium, after which the 
residue undergoes statistical decay. 
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Fig. 9. Lower frames: distributions of observed total thermalized energy per event for 
^He + ^^^Au(left) and ^He + "''^Ag (right); upper frames: correlation between total 
thermalized energy and transverse energy. Error bars indicate standard deviations 
of distribution widths (±o") and are representative of data. Systems are defined on 
figure. 

Further, the BUU calculation that incorporates clusters, a momentum-dependent 
potential, in-medium cross sections and a time scale of t ~ 65 fm/c is able to 
describe the excitation-energy probability distribution and the average source 
Z and A as a function of E*/A. One important aspect of this code is that the 
inclusion of clusters serves to enhance energy deposition in the heavy residue. 
This feature is missing in other codes designed for transport calculations with 
GeV hadron beams. 



3.3 Sideways Peaking 



Finally, the ISiS data demonstrate the important role that reaction dynam- 
ics exert on the statistical decay properties of hot residues. Earlier inclusive 
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Fig. 10. Normalized yields of nonequilibrium light charged particles as a function 
of the source excitation energy for the p+Au reaction at 14.6 GeV/c. Circles repre- 
sent data and lines represent filtered ratios from different versions of the transport 
calculations, at different times in the reaction. 

studies [42,43,44] showed the existence of sideways peaking of IMF angular 
distributions for p + A reactions above ~ 10 GeV. This result was interpreted 
as possible evidence for dynamical IMF emission during the initial phases of 
the cascade, possibly signaling the existence of nuclear shock wave effects. Ex- 
clusive IMF angular distributions obtained with ISiS support a more mundane 
origin for the sideways peaking. 

In [45] it is confirmed that the sideways peaking develops only above beam mo- 
menta of 8-10 GeV/c and that degree of peaking increases as IMF multiplicity 
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Fig. 11. Dependence of relative angular distributions on IMF kinetic energy for Z = 
5-9 fragments formed in events with IMF multiplicity M > 4 for 5.0 GeV/c (left) 
and 14.6 GeV/c proton (right) beams incident on ^^'^Au [45]. Angular distributions 
are normahzed to 1.0 at 160°. IMF kinetic-energy bins are indicated in the figure. 
Solid hues in the left and center panels are results of moving-source fits. Right-hand 
panel shows prediction of an INC/SMM calculation [31,46] for the 14 GeV/c p + 
"'^^'^Au reaction, binned the same as the data. 

and charge increase. This result indicates that the peaking is associated with 
high deposition energy collisions, where multifragmentation is the major decay 
mode. In Fig. 11 relative angular distributions for Z = 5-9 fragments produced 
in yiiMF > 4 events are compared for reactions with 5.0 GeV/c vr" (left panel) 
and 14.6 GeV/c p (center panel) beams on a ^^^Au target. Fragment kinetic- 
energy cuts of E/A = 1.2-3.0, 3.0-5.0 and 5.0-10.0 MeV are imposed on the 
spectra and all angular distributions are normalized to unity at 160°. (In this 
regard it should be stressed that most of the cross section is concentrated in 
the kinetic energy bins below 5 MeV.) It is observed that as the IMF kinetic 
energy decreases, the angular distributions become increasingly isotropic for 
both energies. However, whereas for the tt" beam the angular distributions 
remain forward-peaked, for the 14.6 GeV/c proton case sideways peaking is 
observed - and the maximum differential cross section evolves to larger angles 
as the IMF velocity decreases. 

The observation that sideways emission is favored by high beam momentum, 
high IMF charge and low IMF kinetic energy suggests a possible origin in the 
kinematics of the residue rather than in dynamical emission. The diffractive 
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nature of the initial N-N collision at GeV momenta preferentially produces a 
secondary nucleon or N* that recoils 70° - 90° to the beam axis, with the angle 
growing as the beam momentum increases. Subsequent dissipation during the 
cascade imparts a transverse velocity component to the heavy residue. The net 
result is that statistical fragment emission from the residue is focused non-axial 
to the beam direction - which affects the lowest energy IMFs most strongly. 
This conjecture is reinforced by INC/SMM (Statistical Multifragmentation 
Model) calculations [31,46], shown in the right-hand frame of Fig. 11 and is 
also illustrated in Fig. 4. Thus, these results, coupled with IMF-IMF angular 
correlations and sphericity and coplanarity distributions, described in [47], do 
not support arguments for dynamical effects such as shock waves as a source 
of sideways-peaking observed in inclusive angular distribution studies. 



4 Statistical Decay: Multifragmentation 

4-1 Calorimetry 



In any attempt to describe a system in terms of thermodynamics, a knowledge 
of the heat content is fundamental. For hot nuclei, this energetic factor is 
expressed in terms of the excitation energy per residue nucleon, E*/A. In this 
section we examine the procedures for determining E* and A for the ISiS data 

[40,48,49]. 

For each reconstructed event, the excitation energy of the emitting source is 
calculated as follows: 

Mcp 
Kource = E Kcp{t) + M^ < K^ > +E, - Q. (4.1) 



Here K^p is the kinetic energy for all thermal charged particles, M„ is the 
multiplicity of thermal neutrons with average kinetic energy < Kn >, E^ 
is the total energy emitted by gammas, and -Q is the removal energy (the 
negative of the reaction Q-value). Each of these terms requires assumptions, 
as described in the following. 

In calculating the charged particle contribution to E* in Eq.(l), cascade/preequilibrium 
emissions prior to thermalization must be removed from the sum. Separation 
on an event-by-event basis is not fully transparent due to the time evolution 
of the cooling process. Fig. 12 presents angle- integrated spectra for Z = 1,2,3 
and 6 nuclei from the 8.0 GeV/c vr" -|- ^^^Au reaction. Because source veloc- 
ities are low (~ 0.01 c), kinematic effects are small. The principal features 
of the spectra are a Maxwellian low-energy component, which we attribute 



to thermal events, and an exponentially-decreasing high-energy tail due to 
nonequilibrium processes. 

The spectra have been decomposed [41,48,49] with a two-component moving- 
source model [50]. In Fig. 12 the thermal source [51,52] is described by the 
dashed lines, the nonequilibrium source [50] by dotted lines and the total 
by the solid line. The nonequilibrium component is most important for the 
hydrogen isotopes, but is also a significant fraction of the He yield. In the 
proton case the two-component model is insufficient to account for the high- 
energy portion of the spectrum, suggesting evidence for three components: 
thermal at low energies, preequilibrium at intermediate energies, and fast cas- 
cade particles at the highest energies. Such a picture is consistent with the 
BUU transport calculations in Sec. 3. As seen in Fig. 13, nonequilibrium emis- 
sion persists, even at the most backward angles. Fig. 13 also illustrates the 
weak-angular dependence of the spectra. For IMFs, the nonequilibrium com- 
ponent decreases in yield as the IMF charge inreases. For carbon and heavier 
IMFs, the angle-integrated preequilibrium yield is negligible. 

Since H and He isotopes constitute most of the charged-particle yield, calorime- 
try requires a systematic procedure for distinguishing between thermal and 
nonequilibrium emissions on an event-by-event basis. This need was a primary 
motivation for the moving-source fits to the spectra. From moving-source anal- 
yses of all the spectra from the 1.8-4.8 GeV ^He -|- ^^^Au reaction [41], it was 
concluded that a sharp cutoff approximation gave a satisfactory account of 
the fraction of thermal events in the spectra. The thermal cutoff values were: 

Kcp{Z = 1) < 30MeV and (4.2) 

KcpiZ > 2) < (9Z + AO)MeV. (4.3) 

These cutoffs are indicated by the arrows in Fig. 12, and correspond approx- 
imately to the break in the slope of the spectral tails. This method of deter- 
mining the thermal yield was also compared with the integrated yield from 
the two-component fits, which led to slightly lower E* values [48]. For the 
calculation of excitation energy in this work, the cutoff values of Eqs. (4.2) 
and (4.3) were employed. The resulting thermal yields are isotropic in the 
center-of-mass frame [48]. 

In contrast the E* analysis of the EOS 1 GeV -'^^'^Au -|- ""^^C data [53] used a 
sharp cutoff assumption of E/A < 30 MeV for all fragments. This approach 
yields significantly higher E* values, as shown in Fig. 14, largely due to the 
inclusion of preequilibrium He ions in the sum of Eq. (1). When this difference 
is taken into account, the ISiS and EOS experiments are in good agreement 
in those areas where they overlap. While the EOS calibration may lead to 
high E* values, the ISiS sharp cutoff approximation may underestimate E* 



19 



Z=l 




e 
s 




- -■ 




\A 








-3 
10 


• 


1 1 


1 1 1 .1 1 1 


II 


- 


C&D 




Z=3 




-2 
10 


-I 




\ V 








- 




\ 


R 


o 




-3 




1 , ,■ 


\ 




l\ll 1 




20 40 60 SO 100 

Kinetic Energy (MeV) 



Fig. 12. Angle-integrated kinetic-energy spectra in the laboratory frame for Z = 1, 2, 
3 and 6, as indicated in each panel, emitted in the 8.0 GeV/c vr" -|- ^^^Au reaction. 
Open points correspond to data. Dashed (dotted) lines represent the thermal-like 
(nonequilibrated) component of the moving-source fit. The solid line is the sum of 
the two fits. Upper cutoff energies [41] are shown by vertical arrows. 
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Fig. 13. Angle-dependent spectra for Z = 1. Reaction is 8.0 GeV/c vr + Au. 

for high excitation energies. This problem is illustrated in Fig. 15 where the 
LCP kinetic energy spectra are plotted for several E*/A bins. The sharp cutoff 
distinction is clear at low excitation energies, but the two components blend 
into one as E*/A increases, blurring the separation. 

The second major uncertainty in determining E*/A for the ISiS data is the 
unmeasured neutron contribution to both the thermal sum in Eq. (1) and 
the cascade/preequilibrium multiplicity as it affects the source mass [40,48]. 
In order to estimate the thermal-like neutron component, we have normal- 
ized the neutron charged-particle correlations reported by [39] to the ISiS 
charged-particle results. The measured correlations, shown in Fig. 16, are rea- 
sonably well described for charged-particle multiplicity M > 4 by model sim- 
ulations [46,54] and show the same qualitative behavior as has been observed 
in heavy- ion reactions [55]. A mass-balance procedure [56] does not work well 
for the ISiS data. The rapid rise in neutron multiplicity at low energies, where 
charged-particle multiplicities are low, makes the ISiS E* values increasingly 
uncertain below E* < 200 MeV. 

Several assumptions have been employed in order to estimate the average 
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Fig. 14. Excitation-energy distributions for sharp cutoff assumptions of Eqs.(4.2) 
and (4.3) compared with a cutoff value of Kcp < 30A MeV. 

neutron kinetic energy as a function of E*/A. For the 4.8 GeV ^He reactions < 
Kn > was estimated from Coulomb-corrected proton spectra and then iterated 
to obtain a consistent value < Kn >= 2Tth, where T^h = {E * /cbY^'^ and a = 
A/11 MeV^^ [40]. For the hadron-induced reactions several relationships were 
explored, including Fermi-gas and It should be noted that the IMF yield is 
actually largest for the E/A = 1.2-3.0 MeV bin and smallest for the 5.0-10.0 
MeV bin. Maxwell-gas assumptions with level density parameters a = A/8 
MeV~^ and a = A/13 MeV""*^ [48]. Comparisons were also made with SMM 
[46] and SIMON [54] evaporation codes. Based on this analysis, the SMM 
predictions were used as a conservative estimate of the neutron kinetic energy 
contribution to E*. Eq.(4.1) is then iterated to obtain self-consistency. This 
procedure produced a somewhat lower neutron kinetic energy contribution 
than in [40]. For both the neutron multiplicities and kinetic energies, the 
use of averages leads to loss of fluctuation information in the final excitation 
energies. 

The sharp-cutoff and neutron assumptions, along with the detector geometry, 
are included in the detector filter. The minor contribution to Eq. (4.1) from 
gamma emission is assumed to be E7 = 2(Mcp + M„) MeV in [40] and E7 
= M(Z > 3) MeV [48]. In order to calculate the removal energy (-Q), the 
charge and mass of the source must be reconstructed. The source charge is 
determined from 



Mneq 
Z source = -^tgi ~ /^ Zi[neq), 



(4.4) 



22 



10 -.-. 



a 
s 

w 

D 



Z=l 



E*/A 



2-4 



4-6 



'' 6-9 



Z=2 



''. °°°o„ 



20 40 60 80 100 

Kinetic Energy (MeV) 



Fig. 15. Angle-integrated kinetic-energy spectra in the laboratory for Z=l and Z=2 
as indicated in each panel, and for three bins in E*/A calculated with the cut- 
off assumptions of Eqs. (4.2) and (4.3) [41]. The symbols correspond to bins of 
E*/A=2-4 (open circles), E*/A=4-6 (fihed triangles), and E*/A=6-9 MeV (open 
triangles). Reaction is 8.0 GeV/c vr" + -^^^Au. 



and the mass from 



Mneq 



A, 



^tgt - Yl Zi{neq)- < Mn{neq) >, 



(4.5) 



where Ztgt and Atgt are the charge and mass of the target, Zj(neq) is the 
charge of the measured cascade/preequihbrium particles, and the nonequilib- 
rium neutron multiplicity is related to the nonequilibrium proton multiplicity 
by < M^""«) >= 1.93M^""9). The assumption for < M^""«) > is consistent 
with BUU calculations [35] and experimental results [59] , and is intermediate 
between the A/Z of the target and experimental systematics [60]. The un- 
measured IMF mass is based on the isotopically-resolved data of [61]. In the 
calculation of E*/A the resultant Q values and neutron multiphcities partially 
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Fig. 16. Relation between the mean neutron multiplicity and the charged-particle 
multiplicity. Solid line corresponds to data points reported for LEAR data by Ref. 
[39]; dotted line shows the SMM calculation, and the SIMON evaporation result is 
given by the dashed line. The open squares give the mean multiplicities estimated 
with the mass conservation assumption [56,57]. 

offset one another in Eq. (4.1); e.g., if the neutron term is over-estimated, then 
the removal energy is reduced, and vice versa. 

Several second-order corrections have been investigated and found to have no 
significant effects, among them: the source velocity (~ 0.01c), source emission 
angle [58], and detector threshold effects. Since ISiS does not measure heavy 
residues, the additional assumption is made that all missing mass and charge 
are contained in a single residue. The residue mass distribution obtained in this 
way is in good agreement with measured results from [53]. The consistency 
of the ISiS calorimetry filter has been tested with SMM [46] and SIMON 
[54] calculations that use the measured source mass, charge and excitation 
energy as inputs to the codes [48]. One final consideration is the effect of the 
exponential decrease in the E*/A probability with increasing E*/A, which 
serves to decrease the excitation energy relative to the bin average. This effect 
is most significant for the highest excitation energy bins when the measured 
distribution is deconvoluted. 

The relative fractions of the excitation energy for the LCP, IMF and neutron 
kinetic energy contributions to Eq. (4.1) are plotted in the top frame of Fig. 
17. Over the full range of E*/A the LCP fraction ranges from 25-30%, and 
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Fig. 17. Relative share of excitation energy for various components of the recon- 
struction procedure of Eq.(4.1) as a function of E*/A for the 8.0 GeV vr" + ^^^Au 
reaction. Top frame: hght-charged particle kinetic energy (solid squares), neutron ki- 
netic energy (open circles), and IMF kinetic energy (open triangles). Bottom frame: 
total particle kinetic energy (open circles) and Q values (solid squares). 



that for neutrons from 20-25%, accounting for over half of the total E*. The 
IMF kinetic energy fraction is relatively small, ranging from negligible values 
at low excitation energies to a near-constant value of ~ 10-12% above E*/A 
^ 6 MeV. In the bottom frame of Fig. 17 the excitation energy fraction due 
to removal energy (-Q) is compared with the fraction for total kinetic energy. 
Except for low E*/A, where the calorimetry is most uncertain, the total kinetic 
energy sum is a near-constant factor of two greater than the removal energy. 

Overall we estimate that the values chosen for the ISiS E*/A data and as- 
sociated assumptions could reasonably be lowered by ~ 5% or increased by 
~ 15%. In the following sections we examine the data as a function of the 
calorimetric E*/A procedure described in this section. 



25 



4-2 Thermal Observables 



Of the many experimental signals for thermal behavior in hot nuclei, the 
most transparent are found in the evolution of the spectra as a function of 
excitation energy. Two fundamental tests that a statistically-decaying system 
must face are: (1) is particle emission isotropic? and (2) are the kinetic energy 
spectra Maxwellian? In Fig. 18 invariant cross sections (parallel vy versus 
perpendicular v_l) velocity components are shown as a function of excitation- 
energy bins for hydrogen and carbon ions measured for the 8.0 GeV/c 7r~ + 
^^'''Au reaction. For energetic hydrogen ions (v > 0.2c) at all E*/A one observes 
a spray of forward-emitted particles that originates in cascade/preequilibrium 
processes. The low-energy part of these plots is nearly symmetric about zero 
velocity, indicating emission from a randomized source moving with an average 
velocity of ~ 0.01c. The ISiS acceptance for thermal-like particles is defined by 
the dashed line in Fig. 18, i.e., the sharp cutoff assumption of [41]. The isotropy 
of the projected angular distributions for thermal LCP and IMF emission, as 
well as the relative insensivity to excitation energy, has been demonstrated in 
[47,48]. In addition, an event-shape analysis for the thermal IMFs from the 
4.8 GeV '^He -|- ^^'''Au reaction is found to be consistent with the existence of 
a randomized system that disintegrates on a very short time scale [47]. 



4-2.1 Fragment Spectra 

The Maxwellian character of the LCP and IMF kinetic-energy spectra is il- 
lustrated in Figs. 12, 13, and 15. Figures 12, 13 and 18 show the systematic 
Coulomb shift in the spectral peaks due to kinematic behavior and fragment 
charge. Figure 12 reveals two opposing trends that become apparent when 
the LCP spectra are gated on E*/A. First, the thermal slopes become flat- 
ter as E*/A increases, the expected result of the increasing temperatures. In 
contrast, instead of showing the expected temperature-dependent increase in 
the spectral (Coulomb) peak energies for a system at normal nuclear density, 
the spectral peak energies decrease. The net effect of these opposing effects 
is that the average mean kinetic energy for thermal particles is essentially 
independent of excitation energy, as shown in Fig. 19. This figure also shows 
the expected increase in the average kinetic energy as a function of fragment 
charge. Within the context of Figs. 15 and 19, the observed fragment mean 
kinetic energies can be interpreted in terms of a compensation between two 
competing factors: an increase in temperature offset by a decrease in the source 
density as the excitation energy increases. This behavior is explored in greater 
detail in Sec. 4.3. 
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Fig. 18. Contour plot of longitudinal vm vs. transverse v^ velocity of hydrogen (left) 
and carbon (right) fragments from the 8.0-MeV/c vr" + ^^"^Au reaction for several 
bins in E*/A. Solid lines indicate geometrical acceptance of the ISiS array; dashed 
line gives the thermal cutoff velocity [48] , not corrected for source velocity. 

4.2.2 Multiplicities 

An important signature of multifragmentation, and its possible relation to a 
nuclear liquid-gas phase transition [3,62,63], is the multiplicity of IMFs in an 
event. The models predict that above E*/A ~ 4-5 MeV, multiple IMF emission 
should appear. Figure 20 examines this prediction, presenting the average IMF 
multiplicity (top frame), emission probabilities for a fixed multiplicity (middle 
frame) and the probability for emitting three or more IMFs relative to two 
or less (bottom frame). The unmeasured heavy residue is not included in 
these probabilities and M is derived from a Monte Carlo reconstruction of the 
measured fragment multiplicity N that accounts for detector geometry and 
thresholds. 

In the top frame the IMF multiplicity is shown to increase monotonically, 
with no apparent deviation near E*/A k. 5 MeV. However, when the averages 
are decomposed into specific probabilities (middle frame), it is observed that 
as E*/A increases, the probability for increasing N/mf opens up systemati- 
cally with increasing excitation energy. The probability for emitting three or 
more IMFs (the classical definition of multifragmentation [3]) is seen to grow 
rapidly near E*/A ~ 4 MeV, so that above E*/A ~ 5 MeV multifragmenta- 
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Fig. 19. Fragment mean kinetic energy as a function of IMF charge calculated in 
the source frame for four bins of excitation energy, as indicated on the figure. Data 
are from the 8.0-GeV/c vr^-induced reaction on ^^ Au. 

tion becomes the dominant decay channel. This growth in IMF multiplicity is 
accompanied by a corresponding growth in the thermal LCP multiplicity. 



4-2.3 Charge Distributions 

Another important aspect of the multifragmentation mechanism is the distri- 
bution of fragment sizes (Z), of relevance to the question of critical phenomena 
and the liquid-gas phase transition. The ISiS charge distributions have been 
analyzed in terms of a power-law function, (t{Z) oc Z~'^ , shown in Fig. 21. The 
results are nearly identical for all of the hadron-induced reactions and behave 
similarly for the ^He -|- ^^'^Au data [47]. At the lowest excitation energies, the 
large values of the power-law exponent r imply that small fragments dominate 
the charge distribution, consistent with lower-energy proton-induced reactions 
[61,64,65]. As the system is heated, r values decrease, signifying the increas- 
ing tendency to form larger clusters. A minimum in r is reached near E*/A ^ 
5-6 MeV, corresponding to the rapid increase in IMF multiplicities. The ten- 
dency to form lighter clusters at high excitation energy is most likely due to 
the dissolution of the larger clusters in the heat bath and/or the formation of 
highly-excited clusters that undergo secondary decay. The average variances 
of the Z distributions have also been measured [66] and the relation of both r 



28 




E /A (MeV) 

Fig. 20. Top: average number of IMFs for the observed yield (cfosed circles), the yield 
corrected for geometry (solid triangles), and the yield corrected for both geometry 
and fragment energy thresholds (open triangles) as a function of E*/A for the 8 
GeV/c 7r~ + -'^^'''Au reaction. Middle: probability for given number of detected IMFs 
as a function of E*/A. Bottom: probability for IMF multiplicity M > 3 (circles) and 
M<3 (triangles). 

and the average variance is discussed in Sec. 6 with regard to phase transition 
arguments. 



4-2.4 Cross Sections 

Comparison of the excitation-energy distribution in Figs. 7 and 14 and the IMF 
multipUcity probabihties in Fig. 20 demonstrates that multifragmentation is 
a small fraction (< 10%) of the total cross section in these light-ion reactions, 
most likely originating in the lowest impact-parameter collisions. For the ^He 
beams the calculated total reaction cross section is approximately 2000 mb 
for the ^^^Au target and 1400 mb for ""^Ag [67,68]. For the hadr on- induced 
reactions, the reaction cross section is about 1900 mb [67,68]. The multifrag- 
mentation {MjMF > 3) cross section ranges from about 50 mb for the 5 GeV/c 
7r~ beam to 100-125 mb at the higher energies, where excitation-energy satu- 
ration occurs. 
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Fig. 21. Top: average ratio of observed and geometry-corrected IMFs per residue 
nucleon as a function of E*/A; symbols are defined in the bottom frame. Bottom: 
power law parameters r from fits to the charge distributions as a function of E*/A 
of the residue. 

In Table 2 the cross section dependence on target, beam energy and IMF mul- 
tiplicity is presented for the ^He-induced reactions. The yields clearly increase 
with target mass and projectile energy, with cross sections that range from 9 
mb for the 1.8 GeV ^He + "■'^^Ag case to 190 mb for the ^^'''Au target at 4.8 
GeV. The effect of excitation-energy saturation in the ""^Ag system near 3.5 
GeV bombarding energy is evident. 



4-2.5 Source Charge 

Determination of the Z and A of the emitting source, as well as the unmea- 
sured heavy residue(s), is also an important component of the reconstruction 
process described in Sec. 4.1. In the top frame of Fig. 22 the average fraction of 
the source charge relative to the target charge is shown as a function of E*/A. 
As the excitation energy increases, the effect of nonthermal particle emission 
becomes quite strong, leading to average source charges of Z ~ 60 at E*/A ~ 
8 MeV. The middle frame of Fig. 22 shows the average missing charge in our 
reconstruction procedure, assumed here to be a single fragment. Comparison 
with the SMM model [46] gives good agreement and these values closely corre- 
spond to the largest fragment distribution observed in the EOS 1 GeV ^^^Au 
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Table 2 

Cross sections for the ^He-induced reactions 



^He + ""*Ag 



^He + ^'^'^Au 



Efteam (GcV) 


1.8 


3.6 


4.8 


1.8 


4.8 


MiMF 






cross section (mb) 






1 


140 


160 


190 


270 


300 


2 


28 


77 


98 


66 


170 


3 


7.4 


26 


28 


14 


110 


4 


1.3 


6.4 


6.5 


2.9 


54 


5 


0.2 


1.2 


1.2 


0.5 


20 


6 


0.03 


0.2 


0.2 


0.1 


6.6 


7 


- 


- 


- 


- 


1.8 


8 


- 


- 


- 


- 


0.4 


9 


- 


- 


- 


- 


0.1 


10 


- 


- 


- 


- 


0.02 


(t{Mjmf > 3) 


8.9 


34 


36 


18 


190 


(Timf{M > 1) 


170 


270 


320 


350 


660 


0-7MF (total) 


220 


430 


490 


460 


1300 



+ C reaction (when corrected for differences in thermal excitation energy). 

[53]. 

Beyond E*/A > 6 MeV the assumed fragment corresponding to the missing 
charge (Zmaxi) is an IMF (Z < 20). Finally, in the bottom frame of Fig. 22 
the average charges of the second {Zmax2) and third {Z^axs) largest fragments 
are shown as a function of excitation energy. The sizes of the second and third 
largest fragments remain nearly constant above E*/A ~ 5-6 MeV, in line with 
the results of Fig. 21, and are also in relative accord with SMM predictions up 
to E*/A ~ 7-8 MeV. At higher excitation energies the data and SMM diverge, 
most likely due to the storage of excess excitation energy in fragments in the 
model, leading to secondary decay. This divergence suggests that the fragments 
are emitted relatively cold, as argued in [69]. 

In summary, the thermal observables from the ISiS data present a picture of a 
system that decays isotropically from a source with velocity ~ 0.01 c. Above 
E*/A ~ 5 MeV multifragmentation {Mjmf > 3) becomes the dominant decay 
mode and the spectra suggest emission from a dilute/expanded source. Near 
this energy the largest fragments are formed. Thus, these signals are qualita- 
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Fig. 22. Dependence of fractional source charge and IMF charges as a function of 
E*/A for the 8 GeV/c vr^ + ^^^Au reaction. Top: fractional source charge. Middle: 
missing charge in ISiS, assumed to be the largest fragment, and SMM prediction 
for missing charge (solid line) and for largest fragment (dashed line), both passed 
through the ISiS filter. Bottom: charge of two largest observed fragments, solid line 
is the SMM prediction for second largest fragment {Zmax2) and dashed line for third 
largest fragment (Zmaxs)- 

tively consistent with expected observables from a liquid-gas phase transition 
and in the following sections we examine this question in greater depth. 



4-3 Breakup Density and Expansion 



A knowledge of the dependence of nuclear density on thermal excitation is of 
central importance to our understanding of nuclear compressibility and the 
equation-of-state of finite nuclear matter. In addition, the breakup density is 
particularly relevant to models of multifragmentation phenomena, which as- 
sume that at sufficiently high temperatures, thermal pressure and Coulomb 
forces drive nuclear expansion and subsequent disintegration of the system 
[3,46,54,62,63]. Perhaps the most direct experimental signal of the breakup 
density is provided by the centroids of the peaks of the IMF kinetic en- 
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Fig. 23. Energy spectra at 43° and 137° for carbon fragments emitted from the 
4.8 GeV '^He + ^^'^Au reaction as a function of excitation energy. Symbols are as 
follows:<E*/A> = 3.4 MeV (•); 4.6 MeV (□); 5.7 MeV (A); 6.8 MeV (O); 7.9 MeV 
(V). 

ergy spectra as a function of E*/A, as mentioned in previous sections and 
in [10,13,47,48,64,70]. 

Breakup densities have been derived from the systematic Coulomb shifts of 
the spectral peaks for a series of IMF data sets that span the excitation-energy 
range E*/A = 0.9 to 7.9 MeV. Inclusive data were analyzed for the 200-MeV 
^He + 19^ Au [71] and E/A = 20-100 MeV ^^N + i97Au [70] reactions. Exclusive 
data were based on the 4.8 GeV ^He + ^^^Au reaction [72]. The spectra for each 
system were measured with low kinetic-energy thresholds and covered nearly 
the entire 180 degree angular range. The threshold/angle criteria are essential 
for obtaining stable moving-source fits, which require very good definition of 
the low-energy component of the spectra. Details of the analysis are discussed 
in greater detail in [72] and [73] . Representative IMF kinetic-energy spectra for 
carbon fragments from the 4.8 GeV ^He + ^^^Au system are shown in Fig. 23 
as a function of E*/A. The decrease in the peak centroids with increasing 
excitation, opposite of expectations for a thermal source at normal density, is 
evident. 

The spectra were analyzed in terms of a two-component (three for ^^N) moving- 
source model [50] consisting of a thermal-source described by a transition-state 
formalism [51,52], a nonequilibrium source that assumes a Maxwellian shape. 
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and for ^^N, a projectile-breakup source. The thermal source, of primary in- 
terest for this analysis, included the following parameters: the source velocity, 
a fractional Coulomb term kc, a spectral slope temperature and a barrier fluc- 
tuation variable. The decreasing Z of the source with E*/A was also taken into 
account. This formalism is designed primarily for binary breakups. For multi- 
fragmentation events, it is assumed that this procedure provides a first-order 
approximation to the Coulomb field that exists between a given IMF and the 
average of the residual nucleons. For density determinations the average frac- 
tional Coulomb parameter <kc> for IMFs is the sensitive parameter and is 
determined relative to fission fragment kinetic energy systematics [72,73]. In 
the top frame of Fig. 24 values of the Coulomb parameter are plotted versus 
E*/A, where kc = 1 corresponds to nuclei at normal density po- Most striking 
about Fig. 24 is the sharp decrease in < fc^ > in the excitation-energy interval 
E*/A ^ 2-4 MeV, suggesting emission from a source with a modified Coulomb 
field. It is also of note that this decrease matches the major changes in other 
reaction variables, as shown in Sec. 4.2. 

In order to obtain the breakup density, it is assumed that for these light-ion- 
induced reactions the thermal source is spherical and expansion is radially 
symmetric. In this case, from Coulomb's law the density expression reduces 
to 

< p > /po = 4, (4-6) 



since < kc > «< re >~^, where yc is the mean separation distance at breakup. 
The result is shown in Fig. 24 where <p >/po is plotted as a function of E*/A. 
Up to E*/A ^ 2 MeV the density appears to correspond to normal density. 
Between E*/A ^ 2 and 5 MeV, Fig. 24 indicates a systematic decrease in 
density from <p >/po ^ 1.0 to <p >/po ^ 0.3. Above E*/A > 4 MeV, a 
nearly constant value of <p >/po ^ 0.30 is found within experimental error. 

In summary, this spectral-shape analysis indicates that above E*/A > 2 MeV, 
nuclear breakup occurs from an increasingly expanded/dilute configuration. 
Beyond E*/A > 4 MeV, a value of <p/po >~ 0.3 is found, consistent with 
the predictions of multifragmentation models [46,62]. Finally, the relative con- 
stancy of <p/po > at high excitation energies suggests that a limiting breakup 
density has been reached [74,75]. 

The energy that drives expansion and subsequent multifragmentation is usu- 
ally attributed to either internal thermal pressure [76] or the response to com- 
pressional forces produced in the early stages of the target-projectile interac- 
tion [77]. As discussed in Sec. 3, model calculations of the reaction dynamics for 
the systems studied in this work provide little or no evidence for compression- 
decompression effects. Instead, the fast cascade creates an initial residue with 
lower than normal nuclear density and high thermal energy. Thus, in highly 



34 




w - 



.^- 



a. 

V 



M, 



1 i 



-i \- 



H 1- 



a 



k 



i i 



-J — 1„ 



EVA CMcVl 



Fig. 24. Top: Dependence of the average fractional Coulomb parameter < kc > 
as a function of excitation energy. Symbols are as follows: 200 MeV ^He (solid 
square);^^N (solid triangle); 4.8 GeV ^He (open square). Bottom: Average density 
<p/po > as a function of E*/A derived from the kc values in the top panel. 

asymmetric collisions at GeV energies, only the thermal pressure and Coulomb 
field are of primary relevance to the subsequent expansion process. 

Two stages of thermal expansion leading to multifragmentation can be schemat- 
ically defined. The first drives the nucleus to the breakup configuration, where 
repulsive Coulomb forces exceed the restoring nuclear force. Although the 
breakup density beyond this point may remain nearly constant, as in Fig. 24, 
the increasing heat content of the source may lead to an additional source 
of radial expansion energy (or flow); i.e., the thermal energy that exceeds the 
minimum necessary to reach the breakup density, defined here as excess expan- 
sion energy, eth- The impact of the excess expansion energy on the spectra will 
be an increased flattening of the high-energy spectral slope, over and beyond 
that expected for the freezeout density and temperature [78]. 

In order to investigate the possible contribution of eth, it is necessary to choose 
a reference point that accounts for thermal and Coulomb contributions to the 
fragment kinetic energies. For this purpose, three statistical models have been 
employed: SMM [46,79], SIMON-explosion and SIMON-evaporation [56]. In 
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Fig. 25. Comparison between experimental and simulated fragment mean kinetic en- 
ergies calculated for two bins in excitation energy. In each panel, data are shown with 
open and solid circles and simulations with dashed and solid lines. The correspond- 
ing bins of excitation energy are indicated on the figure. SMM and SIMON-explosion 
calculations have been performed without additional expansion energy. 

order to minimize any contamination from possible preequilibrium emissions 
in the data, we compare the calculations with the measured average frag- 
ment kinetic energy, for which preequilibrium contributions to the tails of the 
spectra are negligible. The inputs to all three models were the same, using 
the source charge, mass, velocity and excitation-energy distributions for the 
reconstructed data (Sec. 4.1), and then passed through the ISiS filter. 

In Fig. 25 the calculations are compared with mean fragment kinetic energies 
as a function of fragment charge for excitation-energy bins, E*/A = 4-6 and 6-9 
MeV. The evaporative model underpredicts the data substantially, although it 
does give reasonable agreement for E*/A < 3 MeV. Both of the simultaneous 
multifragmentation models describe the mean kinetic energies, as well as the 
multiplicity and charge distributions, for E*/A = 4-6 MeV bin. However, for 
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the E*/A = 6-9 MeV bin both models fall below the data. This shortfall is 
attributed to the existence of excess thermal expansion energy. Using the SMM 
model as a reference point, the excess expansion energy is extracted from the 
difference between the model and the data. In Fig. 26 the results are plotted 
versus E*/A. This analysis indicates that the eth threshold occurs near E*/A = 
4 MeV and then gradually increases to 0.5 A/mf MeV at E*/A = 9 MeV. This 
amount of energy, while small, must be taken into account when performing 
the calorimetry (Sec. 2). In contrast, central heavy-ion reactions show a much 
more dramatic increase in the excess expansion energy, also shown in Fig. 26, 
suggesting that the excess expansion energy observed in heavy-ion collisions 
may be related to the dynamical stage, perhaps due to initial compression. 

To summarize, moving-source analyses of the IMF spectra show a system- 
atic downward shift in the Coulomb peaks, supporting a picture in which 
the breakup density decreases as a function of excitation energy, even after 
correcting for nonequilibrium charge loss. While the density remains nearly 
constant above E*/A = 5 MeV, the flattening of the spectral slopes suggests 
a small, but measurable excess breakup energy that increases nearly linearly 
as eth/AiMF = (0.1 E*/A)-0.4 MeV. 



4-4 Breakup Time Scale 



Central to any interpretation of multifragmentation events in terms of a liquid- 
gas phase transition is the question of time scale. For evaporative cluster 
emission from the liquid phase at low excitation energies, fragments are pro- 
duced from the surface via a binary sequential decay mechanism. This process 
requires relatively long emission times of order 1000 fm/c at low excitation 
energies [80]. In contrast, when the spinodal boundary of the phase diagram 
is crossed, the system falls apart on a near-simultaneous time scale via bulk 
emission from the entire nuclear volume. 

Information about the emission time scale can be extracted by means of the 
intensity-interferometry technique, which probes the mutual Coulomb repul- 
sion between fragment pairs emitted in proximity to one another in space and 
time [81,82,83,84]. This technique constructs the correlation function R for 
fragment pairs of reduced velocity Vred, 

i?(v.ed) + 1 = C ^— ^J-'^^ . (4.7) 

^^ uncorr l^ *^red ) 



^corr is the measured coincidence yield, while '^uncorr is the uncorrected yield 
calculated with the event-mixing technique [83], and the normalization C is 
performed relative to the integral yields of each [85] . 
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Fig. 26. Upper panel: See Fig. 20. Lower panel: comparison between 8 GeV/c 
vr^+^^^Au reactions and central heavy-ion collisions [78]. The shaded area corre- 
sponds to the ISiS excess expansion energies extracted with SMM at 3Vo (upper 
limit) and 2Vo (lower limit). The dashed line summarizes the excess expansion en- 
ergies extracted in central heavy-ion collisions with various assumptions regarding 
the source characteristics. 

The reduced velocity is given by 



Vred 



I Vl - V2 

VZi + Z2 



(4.8) 



where Vj and Zj are the laboratory velocity and charge of the fragments, re- 
spectively. The denominator permits comparison of different IMF Z values. 

Experimental IMF-IMF (4 < Z < 9) correlation functions from hadron-induced 
reactions on ^^^Au between 8.0 and 11.2 GeV/c are shown in Fig. 27 for sev- 
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Fig. 27. Reduced-velocity correlation functions generated for four different excitation 
energy energy per nucleon bins. IMF kinetic energy acceptance in the source frame 
is E/mf/A = 1-10 MeV. 



eral excitation-energy bins. Pairs emitted in close proximity to one another in 
space and time (low ^red) experience a supression in yield due to their mutual 
Coulomb interaction (Coulomb hole). Between E*/A = 2.25 and 5.0 MeV the 
Coulomb hole increases, followed by a nearly constant supression at higher ex- 
citation energies. This effect is in qualitative agreement with heavy-ion studies 
[86,87]. 

In order to extract the emission time scale, an N-body Coulomb trajectory 
calculation [84,88] has been performed, using the experimental source and fi- 
nal product properties as input [85]. The only adjustable parameters in the 
simulation are the source volume, or separation distance between the residue 
and the fragments. The filtered output of the simulation must reproduce both 
the small- and large-angle correlation data, as well as the fragment charge 
distribution and kinetic energy spectra. These conditions impose a significant 
constraint on space-time ambiguities in the simulation. For purposes of calcu- 
lating the Coulomb energy, the separation distance is defined as 
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Ri2 = ro{Ali^ + A'/^f) + d 



where tq = 1.22 fm. Values of d between 2-6 fm provide the best fits to the data 
and are consistent with the density resuhs described in the previous section. 
The emission time t is assigned via an exponential probability distribution, 
g*/*o^ where to is the decay lifetime. 



In Fig. 28 the experimental correlation functions are compared with simula- 
tions for a range of d and time values that yield minimum chi-squared values. 
Between E*/A= 2.0-2.5 and 4.5-5.5 MeV the emission time decreases by an 
order of magnitude, from ~ 500 fm/c to 20-50 fm/c. Above E*/A ~ 5 MeV 
the emission time becomes very short and nearly independent of excitation 
energy, consistent with a near-instantaneous breakup scenario. Similar results 
have been shown for the 4.8 GeV ^He -|- ^^^Au reaction, as shown in Fig. 29 
[89]. Here the experimental gating is performed for all fragments from events 
with E*/A > 5 MeV. Again the results are consistent with breakup times 
between 20-50 fm/c for that fraction of the yield that falls in the multifrag- 
mentation regime. 



The lower panel of Fig. 30 presents the best-fit decay times for events in which 
two or more IMFs are emitted in hadron-induced thermal multifragmentation 
of ^^^Au nuclei. For reference the individual IMF multiplicities are shown in the 
upper panel. The decay lifetimes at low excitation energy are consistent with 
an evaporative mechanism, while at higher energies the very short lifetimes 
support a near-simultaneous breakup. The shaded band in Fig. 30 covers the 
range of space-time values that provide a consistent fit to all of the observables. 
Also shown in Fig. 30 are results for heavy-ion reactions [90], which yield 
similar results, but somewhat longer lifetimes at low energies. 



In summary, the time scales derived from the intensity-interferometry anal- 
ysis demonstrates the evolution from the evaporative to near-simultaneous 
breakup regime. As with the multiplicity, spectra and density evolution dis- 
cussed in Sees. 4.1-4.3, the time scale determinations provide a strong case for 
an interpretation in terms of a transition from surface to bulk emission in the 
excitation energy interval between E*/A ~ 3-5 MeV. 
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Fig. 28. Correlation function for Z = 4-9 IMFs as a function of reduced velocity 
(open circles). IMF kinetic energy acceptance in the source frame is E/mf/A = 
2-10 MeV. Data are gated on source E*/A = 2.0-2.5 MeV (top), 4.5-5.5 MeV 
(center) and 8.5-8.5 MeV (bottom). Solid and dashed lines are results of a Coulomb 
trajectory calculation for fit parameters indicated on the figure. 

5 Thermodynamics 

5.1 The Caloric Curve: Isotope-ratio Temperatures 



One of the most stimulating early results of multifragmentation studies was 
the excitation energy versus temperature curve, or caloric curve, proposed 
by the ALADIN group [91]. By plotting temperatures derived from double- 
isotope ratios [92] as a function of excitation energy, a result was obtained that 
resembles the heating of liquid water to the vaporization phase. Subsequent 
experiments, including those described in this section, produced similar results 
[75] . From the systematic behavior of these data, Natowitz has derived a value 
of the critical temperature of 16 ± 1 MeV for infinite nuclear matter and a 
nuclear compressibility constant K = 232 ± 30 MeV 
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Fig. 29. Reduced- velocity correlations as a function of reduced velocity for the 4.8 
GeV ^He + ^^^Au reaction (points) [89]. Data were selected for pairs of events in 
which N^/j > 11 and kinetic energy (E/A)/mf = 0.7-3.0 MeV and are shown for 
Z = 3,4 fragments (upper frame) and Z = 5-9 fragments (lower frame). Lines are 
results of an N-body simulation with p/ pq = 0.25 and maximum residue size, Z^es 
= 12. Time scales are indicated in the figure. 

In order to construct the heating curve for the ISiS data, the heat content 
was based on the calorimetry described in Sec. 4.1. The double-isotope-ratio 
technique for defining temperature is limited for the ISiS data due to the 
high thresholds for isotope identification. The isotope-ratio temperatures T 
were calculated according to the prescription of Albergo [92] , with correction 
factors K proposed by Tsang [94], 



1/T 



In(aR) — In/t 
B ■ 



(5.1) 



Here B is a binding-energy parameter, a is a statistical factor dependent on 
ground-state spins, and R is the double-isotope ratio. For ISiS, useful isotope 
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Fig. 30. Upper frame: IMF probability as a function of E*/A. Bottom frame: emis- 
sion time as a function of E*/A. Open points are ISiS data under two extremes of 
fitting procedure. Solid points are heavy-ion data from [90]. Parameters for expo- 
nential fits to the data are given on the figure. 

identification is restricted to LCPs so that the relevant ratios are 



Rpd-He = {p/d)/{^He/^He) and 



(5.2) 



Rdt-He = {d/t)/CHe/'He), 



(5.3) 



where all ratios involve only thermal LCPs. 

The definition of thermal LCPs, as discussed in Sec. 4.1, can lead to variabil- 
ity in the value of T that is obtained. For p/d and d/t ratios, there is only a 
small sensitivity to particle kinetic energy [49,72], so the distinction between 
thermal and preequilibrium particles is of minor significance due to the loga- 
rithmic nature of Eq. 5.1. In contrast, as shown in Fig. 31, the '^He/^He ratio 
increases strongly with kinetic energy [40,49,72]. Thus the controlling factor 
in determining T is the ^He/^He ratio, which is dependent on the thermal 
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Fig. 31. He isotope ratios [40] as a function of He energy observed at 43 and 137 
degrees for the 4.8 GeV ^He reaction on Ag (left panel) and Au (right panel). Lines 
are INC/EES model predictions [62,95] for 137 degrees (soUd ) and 43 degrees 
(dashed). Error bars are statistical only. 

cutoffs assumed in the calorimetry. 

In Figs. 32 and 33 the temperature versus E*/A curves are shown for the 4.8 
GeV ^He + ^^^Au system [40] and the 8.0 GeV/c vr" + ^^^Au system [49], 
respectively. For the ^He-induced reaction, which uses only the dt-He ratio, 
there is no plateau, but a slope change is observed above E*/A ^ 2 MeV. 
While the temperature increases from about T = 5-7 MeV in the E*/A ^ 
2-10 MeV range, it deviates markedly from a simple Fermi gas prediction 
(dotted curves in Fig. 32). 

Also shown in Fig. 32 are comparisons with INC/EES (Expanding, Emitting 
Source) [30,62] and INC/SMM [30,46] models. SMM comparisons assume the 
fragments are emitted cold; comparisons with model parameters that produce 
hot fragments deviate strongly from the data [49] . For both model comparisons 
the solid lines are predictions with the experimental cuts defined in Sec. 4.1 
imposed on the model spectra. The results provide fair agreement with the 
data. With the experimental cuts removed (dashed curves), the caloric curves 
are lowered by ~ 1-1.5 MeV per nucleon, yielding better agreement with other 
caloric curves. The difference between the results with and without the exper- 
imental cuts on the model is traced to the fact that the thresholds for isotope 
identificaton in ISiS fall above the spectral peaks, where the ^He/^He ratio is 
much lower [96,61] (i.e., R is larger and T is smaller). Also shown in the right 
frame of Fig. 32 is the thermodynamic temperature predicted by the SMM 
model. 
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Fig. 32. (d/t)/(^He/'^He) isotope-ratio temperature vs. reconstructed E*/A for the 
4.8 GeV ^He + ^^^Au reaction. Left frame compares data with the INC/EES model 
[62,30] and right frame compares with the INC/SMM model [46,30]. Solid curves are 
model predictions with experimental cuts imposed on H and He kinetic energy spec- 
tra. Dashed curves show the effect of removing the experimental cuts. Dotted curves 
show Fermi gas behavior with a=ll MeV~^. For the SMM case the dot-dashed curve 
gives the thermodynamic temperature of the source. 

In Fig. 33 the caloric curves for the tt" + ^^^Au reaction are compared for 
the pd-He and dt-He thermometers. The top frame shows the results without 
the Tsang [94] corrections. With the correction applied, both thermometers 
demonstrate a break in the curve above E*/A ^ 3-4 MeV, with the pd-He 
ratio yielding a more distinct plateau. At the highest excitation energies for 
the dt-He case there is an indication of an upturn in the caloric curve above 
E*/A ^ 10 MeV, suggestive of possible entrance into the vaporization regime. 
Similar results have been observed in other experiments and are also seen in 
the 4.8 GeV ^He data, although the effect occurs at excitations where statistics 
are low and fluctuations large. However, this behavior is absent with the pd- 
He thermometer. Thus, while the high E*/A dt-He results are intriguing, they 
are not convincing. 

The differences in temperature between the slope/plateau regions for the 4.8 
GeV ^He and 8.0 GeV/c n~ reactions can be traced to two factors. First, the 
kinetic-energy thresholds were lower in the former case and second, the energy 
acceptance bins were not quite the same. The net effect is that the ^He/^He 
ratio is lower for the 4.8 GeV ^He measurements; i.e., R is larger and T lower. 
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Fig. 33. Caloric curve for 8 GeV/c ir' + 197Au from (p/d)/ (^He/^He) and (d/t) 
^He/^He) thermometers using measured yields to calculate temperature (top panel) 
and temperatures corrected for secondary decay (bottom panel). 

The temperature dependences on '^He/^He ratio can be used to track the 
evolution of the de-excitation process leading up to thermalization, under the 
assumption that the most energetic emissions are emitted earliest [97,98,99]. 
Such a "cooling curve" is shown in Fig. 34 for the 8.0 GeV/c tt^ + ^^^'^Au 
reaction, where Coulomb-corrected 10-MeV-wide bins have been placed on 
LCP spectra [49]. Note that the higher energy bins correspond to the hard 
exponential tails of the spectra in Fig. 12. The corresponding T vs. E*/A 
calculation reveals a systematic decrease in the isotope-ratio temperatures 
as the kinetic energy bin for the LCPs decreases. The observed sequence of 
caloric curves can be interpreted as evidence for the preequilibrium cooling 
stage between the initial cascade and thermalization stages. An alternative 
explanation is provided by time-dependent EES model, for which particles 
are emitted sequentially from an expanding, cooling source [62]. 
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Fig. 34. The caloric curves for the 8 GeV/c vr^ + ^^^Au reaction from the 
(d/t)/(^He/^He) thermometer, corrected for secondary decay, using four different 
kinetic energy acceptances as given in the graph. The shaded area represents the 
caloric curve region extrapolated to the evaporative region of the H and He spectral 
peaks. 



Due to the systematic trends of Fig. 34 and the lack of mass-resolved data 
in ISiS in the Coulomb-peak region below E/A = 8 MeV, the isotope- ratio 
temperatures would be lower if extrapolated to the thermal LCP peak yields. 
To examine this correction, a linear fit was performed on the cooling curves 
of Fig. 34 and then extrapolated to the Coulomb peak region for element- 
identified LCPs, shown in Fig. 12. The shaded area in Fig. 34 shows this 
extrapolation. Figure 35 compares the caloric curves from similar systems with 
the ISiS data. The left-hand frame shows the difference between the observed 
and corrected ISiS results and emphasizes the sensitivity of the isotope-ratio 
thermometer to the energy acceptance for the LCPs. The center- and right- 
hand frames compare the ISiS results with the ALADIN peripheral Au + Au 
[91] and EOS Au -|- ^^C results [18,53] respectively. Reasonable consistency is 
observed with the corrected ISiS caloric curve and the systematics of Natowitz 
[75]. 
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Fig. 35. Left: Summary of caloric curves from the ISiS data in Figs. 32-34, compared 
with ALADIN data [91] (center) and EOS data [18,53] right. 

5.2 The Caloric Curve: Density- Dependent Fermi- Cas Temperatures 



The density determinations described in Sec. 4.3 suggest an alternative ap- 
proach to measuring the nuclear temperature. Inherent in the Fermi gas model 
is the first-order relationship 



E*/A 



T^ /Px2/3 ^*(Po) 

K{p) po m*{p) ' 



(5.4) 



where K(p) is the density-dependent inverse level density parameter (1/a) 
and m* is the effective mass. This predicted dependence of temperature on 
density provides a method for determining nuclear temperatures, independent 
of isotope ratios. If one assumes that the effective mass ratio is near unity at 
these high excitation energies, then the ratio of K(p) to Kg, the value of K at 
normal density, becomes 



Kip) = Ko{p/poy^'=Ty{E*/A). 



(5.5) 



From fits to data below E*/A < 2.0 MeV, an empirical inverse level density 
parameter of Kq = 11.3 MeV for a density-independent Fermi gas is obtained, 
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Fig. 36. First-order density-dependent Fermi-gas temperatures as a function of ex- 
citation energy. The dashed line is the normal-density Fermi-gas prediction, with T 
= \J\\.?>{E * jA) MeV. Symbols are the same as in Figure 24. 



shown as a dashed curve in Fig. 36. Using this value, Eq. 5.5 becomes 



T{MeV) = [ll.3{p/por'{EyA)] 



2/3^77* //I ^11/2 



(5.6) 



When the average densities derived from the IMF spectra [Sec. 4.3 and Fig. 
24] are inserted into Eq. 5.6, the resultant temperatures produce the caloric 
curve shown in Fig. 36. Up to E*/A ^ 2 MeV the temperature rises according 
to Fermi gas predictions for nuclei at normal density. In the region E*/A = 2-5 
MeV a distinct slope change occurs, corresponding to the decrease in breakup 
density of the emitting source. Above E*/A ^ 5 MeV, Eq. 5.6 with a constant 
value of p/po ~ 0.30 predicts a simple gradual increase in temperatures given 
byT 



2.2 J E* /A. Overall, the density- dependent Fermi gas model yields 
a caloric curve that is strikingly similar to other caloric curve measurements 
for similar reactions [91,53], as well as the corrected caloric curves from ISiS, 
shown in Fig. 34. This result is consistent with statistical model calculations 
that assume p/po ~ 1/3 at breakup [3,46,63] and in qualitative agreement with 
the metastable mononucleus model of Sobotka [100]. A second-order analysis 
[73] of the density data, in which expansion energy is taken into account, 
shows a 1-2 MeV dip in the plateau near E*/A ^ 5 MeV, but is otherwise 
similar. Overall, the temperatures derived from the density-dependent Fermi 
gas model and those from double isotope ratios [94] are in agreement, lending 
additional support to the concept of caloric curve behavior for hot nuclei. 
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Fig. 37. Panel shows Cv and x^ calculated from data. The minimum of x^ and the 
maximum of C^ coincide within experimental uncertainty. 

5.3 Heat Capacity 



In heavy-ion studies evidence for a negative excursion in the heat-capacity ver- 
sus excitation-energy curve has been presented by the MULTICS-MINIBALL 
Collaboration [101]. Based on thermodynamic considerations, this result pro- 
vides possible evidence for a first-order phase transition [102]. The ISiS data 
have been examined for such an effect [103] and the results are shown in Fig. 
37. A sharp negative deviation in the heat capacity Cv is observed near E*/A 
^ 4-5 MeV, consistent with the heavy- ion results [101]. Within experimental 
uncertainties the minimum in the x-squared distribution coincides with the 
maximum in the Cv curve. 

The first-order phase transition argument is reflected in the observations in 
earlier sections that show near E*/A k, 4-5 MeV there is a sharp increase in 
the IMF multiplicity, a rapidly decreasing emission time and density, the onset 
of excess expansion energy, and a distinct slope change in the caloric curve. 
The order of the phase transition is discussed further in the following section. 

In summary, both the corrected isotope-ratio and density-dependent Fermi-gas 
temperature versus excitation energy plots show quasi-caloric curve behavior 
with a transition in the region E*/A ^ 3-5 MeV. Rather than a plateau, the 
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resultant caloric curves show a gradual increase in temperature with added 
heat, but are well below normal Fermi-gas expectations. By gating on LCP 
kinetic energies, it is also possible to construct a " cooling curve" that describes 
the evolution of the reaction mechanism from the cascade step to thermaliza- 
tion. In addition, a signal for a negative excursion in the heat capacity curve 
supports heavy-ion data [102] that are interpreted in terms of a first-order 
phase transition. 



6 The Liquid-Gas Phase Transition: Scaling Law Behavior 



Given the overall agreement of the ISiS data with the qualitative expectations 
for a liquid-gas phase transition, several further questions of a more quan- 
titative nature arise. For example, do the cluster size distributions behave 
according to statistical expectations? How well do the data conform to scaling 
laws for a phase transition? If so, what are the critical parameters and what 
is the order of the phase transition? 

Statistical behavior is an important question, not a priori obvious for systems 
that evolve as rapidly as those formed in GeV hadron-induced reactions. At 
lower energies, where statistical concepts are more appropriate, cluster emis- 
sion probability distributions can be well-described in terms of a binomial 
distribution [104], 

777 I 

P^(E*) = — --tP"(1 - p)'"-". (6.1) 

^^ n\{m-n)\ ^ ^' ^ ' 



Here,n is the IMF multiplicity, m is the number of chances to emit an IMF, 
and p is the binary elementary probability. The values of p and m can be 
extracted from the experimental average multiplicity < N > and its variance, 

<N>=mp and al =< N > {1 - p). (6.2) 



At lower energies it has been found that p is a function of excitation energy, 
giving rise to the concept of thermal scaling [104]. In order to test whether the 
ISiS results follow this same statistical pattern, the 8.0 GeV/c 7r~+^^'^Au data 
have been fit with Eqs. 6.1 and 6.2 [38]. The binomial distribution analysis 
gives very good agreement with the data up to Njmf = 5, as shown in the 
multiplicity distributions in Fig. 38. Some deviations appear for Njmf > 6, 
where statistics become increasingly poor. To investigate the possible depen- 
dence of the binomial parameters on collision violence, the reciprocal of the 
probability factor p was plotted versus the total transverse energy, thermal 
transverse energy and E*/A. The total transverse energy Ef diverges strongly 
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Fig. 38. The experimental (symbols) and calculated (lines) n-fold IMF probabil- 
ity distributions as a function of E*/A. The lines assume a binomial probability 
distribution according to Eq. 6.1 and p and m were extracted from Eq. 6.2. 



Table 3 

Values of the binomial parameter m, the primary source size A^ 

charge Zobs, and the Qvaiue for various E*/A bins. 



the observed 



E*/A (MeV) 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 

m 3.36 4.82 5.61 6.29 6.72 7.67 8.31 7.79 

Zsrc 74.3 71.2 68.4 65.8 63.5 61.1 59.0 57.1 

Aobs 8.75 15.2 21.9 28.1 33.5 38.1 42.2 45.4 

Q„a/«e(MeV) -160 -204 -249 -295 -341 -383 -420 -451 

at high energies, reflecting the contribution of preequilibrium processes to the 
yield. Removal of the preequilibrium component yields an improved scaling 
flt to the data for the thermal transverse energy. The strongest correlation is 
found when scaled as a function of E*/A, where a nearly linear dependence is 
observed, supporting an interpretation in terms of thermal scaling [38]. 

The evolution of the parameter m with excitation energy is shown in Table 3, 
where correlations with the size of the source Zgrc, total observed charge Zobs 
and the removal energy (-Q). The one variable that tracks most closely with m 
is the removal energy, suggesting that m represents an energy constraint that 
allows only certain partitions. 
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The agreement with the thermal scahng concept indicates that p can be ex- 
pressed by a partial decay width 



hjj, 



o 



where uJq is interpreted as the frequency of assault on the barrier B at tem- 
peratures T [38,104]. Defining the intrinsic emission time as tg = l/^o^ the 
emission time is given by 

t = toe^/^,or p = to/t. (6.4) 



Using the emission times derived in Sec. 4.4, a plot as a function of the ther- 
mal scaling variable (E*/A)^/^ shows nearly linear behavior, as seen in the top 
frame of Fig. 39. The bottom frame of Fig. 39 shows the relation between the 
lifetime t and 1/p. A simple linear relationship is observed down to emission 
times of 20 fm/c near E*/A ^ 6 MeV. The evolution of the inverse probability 
1/p at higher energies appears to be independent of time, indicating a mecha- 
nism charge that favors a space-like (bulk) emission scenario, rather than one 
that is sequential, as is the case at lower excitations. 

The apparent statistical nature of the thermal component of the data justifies 
further investigation of expected liquid-gas phase transition properties. One 
method for extracting information relevant to this issue is through a moment 
analysis of the fragment charge distributions [1,105]. Calculations with both 
percolation and statistical multifragmentation models predict that the relative 
variance 72 of the charge distributions will exhibit maxima of 72 ~ 2.1-2.3 near 
the critical point. Brzychczyk [66] analyzed the relative moments of the 4.8 
GeV ^He + ^^^Au reaction and found a variance of 72 = 2.3 ±0.1 near E*/A 
^ 5.5 MeV. Thus, the ISiS data are in good agreement with phase transition 
models and provide further consistency with arguments for a liquid-gas phase 
transition and possible critical behavior in hot nuclei. 

Berkenbusch et al. [106] carried out a global percolation analysis on the 10.2 
GeV/c p + ^^^Au data from ISiS. The bond-breaking probability for the model 
is determined from the excitation energy via the relation [107] 

E,{E*) = 1 - y ^r(3/2, 0, B/T{E*)) (6.5) 



where F is the generalized incomplete gamma function, B is the binding energy 
per nucleon of the source and T is the source temperature, determined from 



the excitation energy T = J E * /a and a = A/ 13 MeV ^ . Input to the model 
utilized experimental values for the excitation energy and source size, with the 
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Fig. 39. Top panel: IMF emission time i as a function of (E*/A)^^'^, from Fig. 
30. The line corresponds to a fit using Eq. 6.4. Bottom panel: plot of 1/p vs time. 
The solid line is a linear fit to the data. The dotted line indicates the "apparent" 
saturation in emission time. 

lattice size fixed by the size of the thermal residue. An important feature of 
the calculation is that account is taken of the secondary decay of the excited 
primary fragments, which is particularly important for the fragile Z = 3-5 
fragments that comprise most of the IMF cross section. 

Figure 40 compares the charge distributions from the data with both filtered 
and unfiltered percolation values. The discontinuity just above Z = 16 is a con- 
sequence of the discrete charge-identification limit in ISiS and the assumption 
that all unmeasured charge resides in a single fragment. Overall, the filtered 
percolation yields are in excellent agreement with the data. 

For values of the control parameter p near the critical value, Pc, the cluster 
number n^ is predicted to scale as 



ns{p) = s ''flip - Pc)S''] (for p ^ Pc), 



(6.6) 



where s is the size of the cluster (Z) and r and a are the two critical exponents 
of percolation theory. The scaling function /has the property that /(O) = 1; 
i.e., the power-law dependence is valid only near p = pc. By associating the 
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Fig. 40. Inclusive charge yield spectra for the reaction p +1^^ Au at 10.2 GeV. The 
round symbols represent the ISiS data. The dotted histogram is the result of the 
corresponding percolation model calculation. The thick histogram represents the 
output of the calculations, filtered through the detector acceptance. 

bond-breaking probability with the temperature and assuming an exponential 
function for /, the fractional IMF yield < n^ > becomes 



<nz>= qoZ^'^expl — ^ ]) 



(6.7) 



where go is a normalization parameter and Tc is the critical temperature. Thus, 
one expects a plot of 



to scale exponentially with (T — Tc)Z^ /T for all fragment sizes. 

The result of a x^ optimization procedure for the theoretical percolation charge 
distributions for this system yields values of the critical parameters a = 0.5 ± 
0.1 and r = 2.18 ± 0.01, in good agreement with percolation analysis for an 
infinite lattice, a = 0.45 and r = 2.18. A similar analysis of the Zjmf = 3-6 
data from ISiS produces the results shown in Fig. 41. The right frame shows 
the inability to achieve scaling when sequential-decay corrections are omitted 
from the analysis. When corrected for sequential decay, as shown in the left 
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Fig. 41. Scaled fragment yields as a function of the scaled control parameter for Z 
= 3, 4, 5 and 6. The left-hand side shows the results of inclusion of secondary decay 
corrections and the right-hand side shows the fit when omitting these corrections. 

frame of Fig. 41, much better scaling behavior is found. The critical exponents 
for this finite system, derived from the data, are a = 0.5 ± 0.1, r = 2.35 ± 0.05 
and (E*/A)c = 5.3 MeV (or Tg = 8.3 MeV for a simple Fermi gas with a = 

A/13 MeV"^). 



Since the infinite size limit of the model contains a continuous phase transition 
for the range of excitation energies covered by the present data set, the scaling 
agreement between the data and percolation theory can be interpreted as 
evidence for a continuous phase transition in nuclear matter. 

The Fisher Droplet Model [108] provides another avenue for examining the cor- 
respondence between IMF emission and a liquid-gas phase transition. Fisher's 
model describes the aggregation of molecules in a vapor into clusters. The 
abundance of a given cluster size A can be written as 



riA = qoA'^exp[{AAfi - CoeA'')/T], 



{6.1 



where in addition to the critical exponents r and a of Eq. 6.7, ua = Na/Aq, the 
number of droplets of mass A normalized to the system size Aq; the difference 
between the actual and liquid chemical potentials is A^ = fi — ^f, Co is the 
zero temperature surface energy coefficient, and e = (Tc—T)/Tc. This equation 
reduces to Eq. 6.7 if A^ = 0. 

One approach taken by Elliott et al. [109] is to modify Eq. 6.8 to take into 
account the Coulomb energy release when a particle moves from the liquid to 
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the vapor phase, which assumes 

UA = qoA-^exp[{AAfi + Ecoui - CoeA'')/T], (6.9) 

where Ecoui is defined as follows 

lA4:{Zsrc — Zjmf)ZimF ., _ye\f a ,r ^ r /r i n\ 

Ecoul = —— — -YTTl X (1 - e '^^)fm - MeV. (6.10) 

ro[{Asrc- AimfY^^ + A/mf] 

Here tq = 1.22 fm and the exponent x^ insures that the Coulomb energy 
disappears at the critical point. Since IMF masses are not measured in ISiS, it 
was assumed that A/Z = 2. Temperatures were determined by the Fermi gas 
approximation of Raduta [110]. IMF acceptance was for IMFs with 5 < Z < 
15, for which preequilibrium effects are small [109]. Subsequent modifications 
of the Fisher model can be found in [111,112]. 

The results of scaling the data according to Eq. (6.9) are shown in Fig. 42, 
where the fragment mass yield is scaled by the power-law prefactor, the chem- 
ical potential term and the Coulomb energy. This quantity is scaled by temper- 
ature and the surface energy parameter A^/T. The scaled data follow Fisher 
scaling over six orders of magnitude, which suggests that this line represents 
the liquid-vapor coexistence line. As supporting evidence, also shown in Fig. 
42 is the scaled cluster distribution from a d=3 Ising model calculation [113] 
for a system that undergoes a phase transition. 

The values of the critical parameters from the ^-squared minimization routine 
are summarized in Table 4 for both the percolation and Fisher scaling analyses. 
Within errors, the Fisher model values of r = 2.28 ± 0.14 and a = 0.54 ± 0.01 
are in good agreement with those from percolation as well as with earlier work 
[114,115]. The surface-energy coefficient Co = 18.3 ± 0.5 MeV is in general 
accord with the liquid drop value of 16.8 MeV. An important result relevant 
to previous scaling analyses is that for the first time it has been possible 
to measure A/i. The measured value of A/x = 0.06± 0.03 substantiates the 
assumption that A/i = in previous analyses. Similarly, the value of x = 1.0 
± .06 insures that the Coulomb energy is small and does not affect the scaling 
significantly. 

Based on this analysis, the phase transition is first order up to the critical 
point of excitation energy E*/A =3.8 ± 0.3 MeV, above which it becomes 
continuous. Analysis of the EOS data [18] yielded a value of E*/A = 4.75 
MeV. However, as discussed in Sec. 3, when the excitation energies are calcu- 
lated with the same assumptions for elimination of preequilibrium particles, 
the ISiS and EOS critical energies are the same. The critical temperature 
for finite nuclei derived from the Fisher scaling analysis is Tc = 6.7± 0.2 
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Fig. 42. The scaled yield distribution versus the scaled temperature for the ISiS 
data (upper) and d = 3 Ising model calculation (lower). For the Ising model, the 
quantity nz/qoA'^'^ /W is plotted against the quantity A'^€/1A35T. Data for T > Tc 
are scaled only as nA/qoA~'^. 



Table 4 

Comparison of percolation and Fisher scaling resutls. 



parameter 


percolation 


Fisher 


r 


2.35 ± 0.05 


2.18 ± 0.14 


a 


0.5 ± 0.1 


0.54 ± 0.01 


T^crit 


8.3 ± 0.2 MeV 


6.7 ± 0.2 MeV 


(E*/A)crit 


5.3 ± 0.3 MeV 


3.8 ± 0.3 MeV 


Co 


- 


18.3 ± 0.5 MeV 


A/x 


- 


0.06 ± 0.03 MeV/A 


X 


- 


1.00 ± 0.06 
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MeV. Based on the Fisher scahng parameters derived from the data, it is then 
possible to construct the two-phase (hquid-gas) coexistence hne over a large 
energy /temperature interval, extending up to the critical point, from which 
the full phase diagram of nuclear matter can be defined [108]. 

In summary, the ISiS data are well- described by binomial reducibility and 
thermal scaling analyses, providing a strong argument for the statistical nature 
of multifragmentation. Further, scaling with the Fisher model can be used to 
define the liquid-gas coexistence line, while a percolation analysis supports a 
continuous phase transition at higher excitation energies. 



7 Summary and Conclusions 



Both the reaction dynamics and the subsequent decay of hot residues formed 
in GeV light-ion-induced reactions on heavy nuclei have been investigated with 
the ISiS detector array. Of primary concern in this effort has been the isolation 
and characterization of multifragmentation events, believed to be the possible 
signature of a nuclear liquid-gas phase transition in finite nuclei. 

Bombardments with proton, antiproton, pion and ^He beams produce an 
exponentially-decreasing distribution of excitation energies that extend up to 
~ 2 GeV in reactions on ^^^Au nuclei. The deposition of excitation energy is 
found to increase as a function of beam energy up to a momentum of about 8 
GeV/c for reactions of hadrons with ^^'''Au and an energy of ~ 4 GeV for the 
^He -|- ""^Ag system. For higher beam energies there is little additional increase 
in deposition energy, presumably due to a tradeoff between beam energy and 
target transparency. Relative to other hadron beams, the optimum projectile 
for achieving high excitation energies is found for 8 GeV/c antiprotons, for 
which the reabsorbtion of some fraction of the decay pions can produce en- 
hanced excitation of the residue. For the same beam momentum, proton and 
pion beams are nearly identical in their excitation-energy distributions. 

The fast LCP component of the spectra has been analyzed with a BUU model 
that incorporates A = 2 and 3 nuclei in the scattering matrix. Best fits to the 
data are obtained with a version that includes modified in-medium scattering 
cross sections and a momentum-dependent potential. As shown in Figs. 4 and 
5, the BUU calculations indicate that for central collisions the hot residues 
are formed in a state of depleted density - due to fast knockout followed 
by preequilibrium processes that occur on a time scale much faster that the 
relaxation time. Simulation of the time evolution of the collision dynamics 
predicts that the entropy per nucleon becomes nearly constant after about 30 
fm/c, suggesting a randomized, but not necessarily thermalized system. From 
comparison of the BUU code with the fast LCP spectra, it is inferred that a 



59 



total time of about 60 fm/c is required to reach a state of quasi-thermalization. 

The most fundamental signature of a liquid-gas phase transition is the obser- 
vation of events in which an equilibrated hot nucleus disintegrates into multple 
LCPs and IMFs. For this purpose, thermal events have been selected that are 
shown to decay isotropically in the center-of-mass system with Maxwellian 
kinetic energy spectra from which nonequilibrium components have been re- 
moved. The cross sections for these events is of order 100 mb. 

Figure 43 summarizes several important features of the data that demonstrate 
a distinct change in reaction mechanism at an excitation energy of E*/A = 4-5 
MeV - all of which support a phase-transition interpretation. In the E*/A = 
4-5 MeV interval the emission of two or more IMFs (accompanied by multiple 
thermal LCPs) becomes the dominant disintegration mode. From a power-law 
fit to the charge distributions, it is found that the fraction of large clusters in 
an average event is a maximum just above the transition excitation energy, 
as predicted by theory. Evidence for the conversion of extra thermal energy 
into enhanced IMF kinetic energy appears near E*/A ~ 5 MeV, although the 
effect is small relative to heavy-ion studies, where compressional heating may 
contribute. 

Two unique results derived from ISiS are the evolution of the breakup den- 
sity and the disintegration time scale as a function of excitation energy. The 
breakup density derived from analyses of the IMF kinetic-energy spectra, 
provides evidence for emission from an expanded/dilute source. The derived 
breakup densities evolve from normal density at low E*/A to a nearly con- 
stant value of p/po ~ 0.3 near E*/A = 5 MeV and above, again consistent with 
the theoretical predictions based on phase transition assumptions. The time 
scale for events with Mjmf > 2 evolves from values characteristic of sequential 
statistical decay at low excitation energy to times of the order of 20-50 fm/c 
for events at E*/A ~ 4 MeV and above, indicative of a near-simultaneous 
decomposition mechanism. 

Similar to previous reports, the data also show evidence for a slope change in 
the dependence of temperature on heat content, suggestive of caloric curve be- 
havior, although a distinct plateau is not observed. Using a density-dependent 
Fermi-gas model to derive temperatures, instead of the conventional isotope- 
ratio thermometer, a caloric curve is obtained that is in good agreement with 
other results. By gating on bins of the preequilibrium spectra, it has also been 
possible to derive a cooling curve for these hot systems. 

Finally, scaling-law fits to the IMF yield distributions provide important con- 
firmation of the statistical nature of the thermal events observed in these stud- 
ies. A Fisher model analysis reveals excellent scaling behavior up to E*/A ~ 4 
MeV, where sequential evaporative emission dominates. This result has been 
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Fig. 43. From top to bottom, the probability distribution for individual IMF mul- 
tiplicities, the power-law exponent that describes the IMF charge distributions, the 
breakup density, excess expansion energy and time scale as a function of excitation 
energy for the reaction of 8.0 GeV/c pions with gold nuclei. 



used to define the two-phase coexistence line for nuclear matter and serves 
as the basis for derivation of the nuclear phase diagram shown in Fig. 44. 
A percolation model analysis that includes IMF secondary decay corrections 
also describes the data well and indicates that at higher excitaton energies, 
the data are described by a continuous phase transition. 

In summary, all of the experimental evidence obtained in these studies is 
consistent with the concept of a nuclear liquid-gas phase transition. However, 
the caveat must be added that given the finite system size and very short time 
scales involved, the classical picture of a phase transition must be invoked with 
caution. It is perhaps better to say that if a nuclear phase transition exists in 
nuclei of a few hundred particles, this is what it looks like. 
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Fig. 44. The reduced density-temperature phase diagram: the thick hue is the calcu- 
lated low density branch of the coexistence curve, the points are selected calculated 
errors, and the thin lines are a fit to the reflection of Guggenheim's equation. 
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